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Using state of the art lattice techniques we investigate the static baryon potential. We employ the multi-hit 
procedure for the time links and a variational approach to determine the ground state with sufficient accuracy 
that, for distances up to ~ 1.2 fm, we can distinguish the Y- and A- Ansatze for the baryonic Wilson area law. 
Our analysis shows that the A-Ansatz is favoured. This result is also supported by the gauge-invariant nucleon 
wave function which we measure for the first time. 



1. Introduction 

Due to the important role which the qq poten- 
tial plays in our understanding of the structure 
of mesons, there exist numerous studies of this 
quantity on the lattice. For a review see ref. 0. 
The three quark potential plays an equally im- 
portant role in the understanding of baryon struc- 
ture. Mass relations between baryons and mesons 
can be made more exact if the form of the bary- 
onic potential is known However, very few 
lattice studies have been made of the baryonic 
potential. Moreover, two such recent studies [|l],D 
have reached different conclusions for the area 
law behaviour of the baryonic Wilson loop: they 
give support to two different Ansatze, called Y- 
and A-law. Since the maximal difference between 
these is a mere 15% for SU{3), a reliable extrac- 
tion of the ground state as well as noise reduc- 
tion techniques are essential in order to resolve 
the dominant area law behaviour Q . 

In this work, employing state of the art lattice 
techniques, we are able to reach sufficient accu- 
racy to distinguish between the two Ansatze up 
to distances of ~ 1.2 fm. At the same time, we 
compare our lattice results directly to the sum 
of qq potentials measured on the same lattices, 
thus avoiding any model assumptions. We also 
compare our procedure and results with those of 
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ref. 1^ in order to understand the different con- 
clusions reached. 

The issue of the dominant area law behaviour 
arises for any gauge group SU{N). Therefore 
it is also interesting to study S'C/(4) and test 
whether the SU(4) results corroborate the conclu- 
sions reached in SU(3). As for SU{3), we choose 
lattice geometries which maximize the difference 
between the two Ansatze which for 5?7(4) is at 
the 20% level. 

As yet a further check, we evaluate the gauge- 
invariant nucleon wave function, using a den- 
sity insertion for each quark line, and examine 
whether it is better described according to the 
Y- or A- law ^. 

2. Baryon Wilson loop 

For the static qq potential the appropriate op- 
erator is the standard Wilson loop. For the static 
baryon potential in SU {N) , the corresponding 
operator is constructed by creating a gauge in- 
variant N quark state at time t = which is an- 
nihilated at a later time T. Explicitly for SU{3), 
the baryon Wilson loop, W^g, shown in Fig. 1, is 
given by 
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P denotes path ordering and r(j) is the path from 
X to y for quark line j. 

The A'^-quark potential is then extracted from 
the long time behaviour of the Wilson loop: 

VW, = - hm ^ln{WN,) . (3) 




Figure 1. The baryonic Wilson loop in SU{3). 
The quarks are located at positions ri, r2 and r^. 

3. Geometries in SU{3) and SU{4) 

Two Ansatze exist in the literature regarding 
the area law behaviour of the baryon Wilson loop: 

• The F-Ansatz: 

In the strong coupling limit, minimization 
of the static energy amounts to giving the 
shortest length, Ly, to the flux tubes join- 
ing the quarks. For SU{3), this is realized 
in general if the three flux tubes meet at 
an interior point known as the Steiner 
point, where their mutual angles are 120''. 
[If one of the angles of the triangle formed 
by the three quarks exceeds 120°, the flux 
tube coming from that summit has length 
zero, and the other two flux tubes meet 
there.] Time evolution of this state in the 
general case produces a three-bladed area 
similar to Fig.|l|, known as the Y- area law. 

For 51/(4) we have more possibilities. 
Minimization of the static energy leads 
to the two stationary solutions shown in 
Fig. ^, namely one configuration with a 




Figure 2. The flux tubes joining four quarks. The 
quarks are located at positions ri,r2,r3 and r4. 
The upper graph shows the local minimum of the 
energy with one Steiner point A, and the lower is 
the minimum with two Steiner points A and B. 

single Steiner point (X-Ansatz), and one 
with two Steiner points A and B {Y- 
Ansatz). The double string between the two 
Steiner points has tension 1.357(29) times 
greater[0 than the other four, single strings. 
If we neglect this difference for simplicity, 
then the F-Ansatz always has lower energy 
than the X-Ansatz. Here, we make no at- 
tempt to distinguish between the Y- and 
X- Ansatze, and we assume that the dou- 
ble string has the same tension as the sin- 
gle strings. Since this assumption has the 
effect of reducing the potential of the Y- 
Ansatz, which is itself lower than in the 
X-Ansatz, it turns out to have no bearing 
on our conclusions. In contrast to SU{3) 
where for any given location of the three 
quarks, the Steiner point and therefore the 
y-Ansatz energy can be computed analyt- 
ically, in SU{4), the two Steiner points in 
the F-Ansatz are obtained by an iterative 
numerical procedure. 

• The A-Ansatz: 

The second possibility for the relevant area 
dependence of the baryonic Wilson loop, 
proposed in ref. ||^, is that it is given by 
the sum of the minimal areas Aij spanning 
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quark lines i and j . Because of its shape in 
SU{3), this Ansatz is known as the A- area 
law. Therefore, we denote by La the total 
length of all interquark distances. 

For SU{3) the maximal difference of 15% be- 
tween the two proposed area laws is obtained 
when the 3 quarks form an equilateral triangle. 

For SU{A) it turns out that this relative dif- 
ference is maximal also for the configuration of 
maximal symmetry among the four quarks. In 
this situation, where the quarks form a regular 
tetrahedron, it reaches 21.96 %. For program- 
ming convenience, we study instead the highly 
symmetric geometries shown in Fig. which all 
give a relative difference of ~ 20%. 
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Figure 3. The geometries considered in S'[/(4). 
For geometry 1 the quarks are placed on a plane 
at equal distances from the origin; for geometry 
2 the positions of the four quarks are along the 
three axes at (1,0,0), (0,1,0) (0,0,1), (0,0,-1); and 
for geometry 3 at (0,0,0), (1,0,0), (0,1,0), (0,0,1). 

We compare our lattice results with the two 
expected forms of the baryonic potential which 
in SUiN) are 
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with Cf = {N^ - 1)/2N and a the string ten- 
sion of the qq potential. Note that in contrast to 
ref. H] we do not allow a to vary. The factor of 
1/(7V-1) in the A- Ansatz makes LaHN-I) < 
Ly always. In addition, we directly compare the 
three- or four-quark potential with the sum of 
two-body potentials measured on the same gauge 
configurations, with no adjustable parameters. 



4. Lattice techniques 

1. Multi-hit procedure: We carried out a com- 
parison between the multi-hit procedure and 
hypercubic-blocking, as proposed in ref. for 
the time links. In brief, hypercubic blocking 
smears the link using staples which all belong to 
the 2^ hypercube surrounding it. Compared with 
the multi-hit procedure, we find that hypercubic 
blocking tends to give larger errors, especially for 
the large Wilson loops. Since accuracy at large 
distances is our objective, we have adopted the 
multi-hit procedure. Moreover, for SU{3) the 
group integral needed to obtain the mean value of 
the link can be computed exactly |l^. We note 
that the reduction factor in the statistical error 
from applying the multi-hit procedure to Wilson 
loops of time extent t grows exponentially with 
t, and that the exponent is further multiplied by 
3/2 in the case of baryonic loops. 

TL. Smearing Correlation matrix: We use a vari- 
ational method to extract the groundstate poten- 
tial. For each quark configuration, we consider 
M different levels of APE smearing, optimized as 
in ref. [Q, and construct an M x M correlation 
matrix C{t) of Wilson loops of time extent t. We 
then solve the generalized eigenvalue problem 



C{t)vk{t) = Xk{t)Cito)vk{t) 



(5) 



taking to /a = 1. We use two different methods 

to extract the groundstate energy: 

1. In the first variant the potential is extracted 

via 



aVk = hmj 



In 
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by fitting to the plateau. Note that with this 
method we also estimate the energy of the first 
excited state. 

2. In the second variant we consider the projected 
Wilson loops 

Wp{t)=v^{to)C{t)voito) (7) 
and fit to the plateau value of —In I Wp{t + 



l)/Wp{t) j . Both procedures gave consistent re- 
sults. The energy of the first excited state was 
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used as a check in the extraction of the ground 
state, ensuring that the contamination is less than 
in the values considered for the plateau. 



the sum of qq potentials extracted from mea- 
surements on the same lattices, i.e. we find 
Vsq w 3/2 Vqq, lu agreement with ref. [||. 



5. Results 

For the baryonic loop in SU{3) we used 220 
configurations at /? = 5.8 and 200 at /3 = 6.0 for 
a lattice of size 16^ x 32 from the NERSC archive. 
For SU (4) we generated 100 configurations at f3 = 
10.9 which gives a similar string tension aa^ as for 
SU{3) at /? = 6.0. 



4.5 r 

4 



V3q 

sum of two-body pots. 



CD 

CD 



cr 
ro r 
> ' 



3.5 



2 

1.5 - 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 
r fm 

Figure 4. The static SU{3) baryonic potential 
at f3 = 5.8 (filled circles). The crosses show the 
sum of the static qq potentials. The curves for 
the A and Y Ansatze are also displayed. The 
quarks are located at {1, 0, 0), (0, 1, 0), (0, 0, 1) and 
r = ri2 = ri3 = r23 = iV^ . 
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Figure 5. As Figure ^ but for /? = 6.0. 

The results at f3 = 5.8 and 6.0 show reason- 
able scaling and are completely consistent with 
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Figure 6. The static SU{3) baryonic potential 
at /? = 5.8 for this work (filled circles) and for 
ref. (crosses) versus the perimeter La of the 
triangle formed by the quarks. 

In Fig. ^ we compare our results with those 
obtained in ref. |l^]. As it can be seen the two 
sets of data are in agreement for small loops which 
correspond to the bulk of the data of ref. . For 
larger loops, the few results of ref. ||l^ tend to lie 
above ours. The statistical errors on these data 
are not quoted but we expect them to be larger 
than ours, especially since the multi-hit procedure 
was not used. The analysis of ref. jl2) differs from 
ours in that we do not allow the string tension to 
vary but take it from the fit to the qq potential. 
So in our approach having fixed the qq potential 
there are no adjustable parameters that enter in 
the two Ansatze. 

Additional support for the A-Ansatz is pro- 
vided by preliminary results on the gauge- 
invariant nucleon wave function, obtained by in- 
serting three density operators at an intermedi- 
ate time t along the three quark lines of a baryon 
propagator. In Fig. |^ we plot the wave function 
versus Ly and La/2. We observe a larger scatter 
when Ly is used as compared to La- The nucleon 
wave function plotted versus La is displayed to- 
gether with a fit to the form cxp(— cL^''^), ex- 
pected if the underlying potential is cx La- As 
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Figure 7. Upper graph: The nucleon wave func- 
tion versus Ly- Lower graph: The nucleon wave 
function versus ^a/S, fitted by a Unear potential 
ansatz. All quantities are in lattice units. 

can be seen, this simple asymptotic form pro- 
vides a remarkably good description of the nu- 
cleon wave function. 

Finally we display the results obtained in 
SU{'l) for geometry 3. They give yet more sup- 
port to the A-area law for the baryon Wilson 
loop. The other two geometries show the same 
behaviour. 

6. Conclusions 

Our results for the static three- and four-quark 
potential in SU{3) and SU{A) are consistent with 
the sum of two-body potentials and inconsistent 
with the Y— Ansatz up to an interquark distance 
of about 0.8 fm. For larger distances, where our 
statistical and systematic errors both become ap- 
preciable, there appears to be a small enhance- 
ment which could be assigned to the admixture 




Figure 8. The static SU{4) baryonic potential in 
lattice units for geometry 3. 

of a many-body component. Nevertheless, for 
the distances up to 1.2 fm that we were able to 
probe in this work, the A-area law gives the clos- 
est description of our data. More refined noise- 
reduction techniques for the large loops will be 
needed in order to clarify or rule out a many- 
body component at larger distances. Preliminary 
results on the nucleon wavefunction also support 
a potential in accord with the A- Ansatz. 

REFERENCES 



G. S. Bah, Phys. Kept. 343 (2001) 1. 
S. Nussinov, Phys. Rev. Lett. 51 (1983) 2081; 
T. D. Imbo, Phys. Lett. B 398 (1997) 374. 
T. T. Takahashi et al, Phys. Rev. Lett. 86 
(2001) 18. 

C. Alexandrou, Ph. de Forcrand and 



A. Tsapalis, |hep-lat/0107006 
C 



Forcrand and 



Alexandrou, Ph. de 

A. Tsapalis, in preparation. 

J. Carlson, J. Kogut and V. R. Pandhari- 
pande, Phys. Rev. D 27 (1983) 233. 

B. Lucini and M. Teper, |hep-lat/0107007 . 
J. M. Cornwall, Phys. Rev. D 54 (1996) 6527. 
A. Hasenfratz and F. Knechtli, 



7. 
8. 
9. 

hep-lat/0103029 . 
10. Ph. de Forcrand and C. Roiesnel, Phys. Lett. 
B 151 (1985) 77. 

M. Guagnelli, R. Sommer and H. Wittig, 
Nucl. Phys. B 535 (1998) 389. 



11 



12. T. T. Takahashi et ai, hcp-lat/0107008 



